Symmetry constraints on phonon dispersion in graphene 
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Taking into account the constraints imposed by the lattice symmetry, we calculate the phonon 
dispersion for graphene with interactions between the first, second, and third nearest neighbors in 
the framework of the Born-von Karman model. Analytical expressions obtained for the dispersion 
of the out-of-plane (bending) modes give the nonzero sound velocity. The dispersion of four in-plane 
modes is determined by coupled equations. Values of the force constants are found in fitting with 
frequencies at critical points and with elastic constants measured on graphite. 
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I. INTRODUCTION 

Since the pioneering experiments on graphene (a single 
atomic layer of graphite) 1,2 , main attention has been de- 
voted to its electronic properties. More recently, Raman 
spectroscopy^ extends to investigations of the lattice dy- 
namics of graphene. It was found that the frequency (sa 
1590 cm -1 ) of the Raman mode in graphene agrees with 
its value in graphite. Also, the overtone of the D mode 
visible almost in all carbon-consisting materials was ob- 
served at about 2600 cm . However, this information is 
very meagre and does not provide a way to describe the 
lattice dynamics. The detailed knowledge of the lattice 
dynamics and electron-phonon interactions^ is needed for 
interpretations of the Raman scattering as well as of the 
transport phenomena. 

Several model s i 6 ' ' 8 ' 9 ' 10 ' have been proposed to pre- 
dict the phonon dispersion in graphene and bulk graphite 
from empirical force-constant calculations. A simplest 
approach assumes the diagonal form of the force-constant 
matrix which contains three constants for the interac- 
tion of an atom with all its nth-nearest neighbor. Thus, 
we meet 12 constants for graphene in the popular 4th- 
nearest neighbor approach or 15 constants in the 5th- 
nearest neighbor onsi£. The number of constants could 
be diminished if the model interactions are use d 13 ' 14 i 15 
or if the phonon dispersion is considered only for the dis- 
tinctive directions in the Brillouin zonei^. 

On the other hand, we can use the most recent 
result a 12 ! 17 ' 18 ! 19 ' 20 ! 21 of the first-principal calculations for 
the phonon dispersion in graphene and graphite. Com- 
parison of that results for the high-frequency modes (see 
Table 1) shows disagreements as large as 50 cm -1 be- 
tween the various approaches. The discrepancies could 
come either from an assumption that the force-constant 
matrix for the atom-neighbor interaction has a diago- 
nal form or from an overestimation of the low-frequency 
modes. It is evident that atoms move more freely in out- 
of-pane direction in graphene than in graphite. There- 
fore, the frequencies of the out-of-plane mode in graphene 
should be less than the corresponding frequencies in 
graphite. Moreover, if the stiffness of the graphene layer 
is neglected, the dispersion of the acoustic out-of-plane 



mode becomes quadratic as seen from the equation of 
elasticity (see also Ref32). The interaction between lay- 
ers in graphite can be estimated from the splitting of 
the low-frequency ZA and ZO' modes in graphite. One 
can see, for instance, from Ref. 2 — that the value of the 
splitting is as much as 130 cm -1 . It means that (i) the 
result of graphene stiffness cannot be larger than that in- 
teraction and (ii) the agreement between the theory for 
graphene and the experimental low-frequency data for 
graphite cannot be better than about 130 cm -1 . 

Here we present an analytical description of the 
phonon dispersion in graphene. This is done within the 
framework of the Born-von Karman model for the hon- 
eycomb graphene lattice including interactions only with 
first, second, and third nearest neighbors and taking the 
constraints imposed by the lattice symmetry into ac- 
count. We show that the out-of-plane (bending) and in- 
plane modes are decoupled from each other. The out-of- 
plane modes are described by three force-constants deter- 
mined in fitting with the Raman frequency and smallest 
elastic constant C44. In the narrow wave- vector interval 
near the T point, the acoustic out-of-plane mode has a 
linear dispersion with the nonzero sound velocity. This 
means that a single graphene layer possesses the small 
but finite stiffness in contradiction with results of Ref. 22 
We should emphasize that the quadratic dispersion of 
the acoustic mode leads to the large contribution (pro- 
portional to the sample size squared of the long-range 
fluctuations, that is much stronger than the logarith- 
mic function in the case of the linear dispersion. Six 
force-constants describing the in-plane modes are found 
in fitting with their frequencies in the critical points and 
elastic constants Cn and C12 of graphite. 



II. PHONON DYNAMICS IN NEAREST 
NEIGHBOR APPROXIMATION 



The equations of motion in the harmonic approxima- 
tion are written in the well-known form 



J2 (»« - a m )< (a m ) 
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where the vectors a n numerate the lattice cells, the su- 
perscripts k, k' note two sublattices A and B, and the 
subscripts i, j = x,y,z take three values corresponding to 
the space coordinates. Since the potential energy is the 
quadratic function of the atomic displacements uf(a n ) 
and uf(sL n ), the force-constant matrix can be taken in 



the symmetric form, $^ B (a n ) = <I>|* A (— a„), and its 



Fourier transform, i.e. the dynamical matrix, is a Her- 
mitian matrix. 

Each atom, for instance, Ao (see Fig. 1) has three first 
neighbors in the other sublattice, i.e. B, with the rela- 
tive vectors Bi = a(l,0), B2,3 = a(— 1, ±V3)/2, where 
a = 1.42 A is the carbon-carbon distance. The sec- 
ond neighbors are in the same sublattice A at distances 
\/3a with the relative vectors Ai^ = ±a(0, VS), A.2,5 = 
±a(-3, a/3) A A 3 , 6 = =Fa(3,\/3)/2. The distance 2a to 
the third neighbors Bi = a(2,0), B 2 3 = a(l,=F\/3) is 
slightly larger. The distance to the fourth neighbors is 
\fla = 2.65a. So, the difference between distances to the 
third and to fourth neighbors is nearly the same as the 
difference between distances to the first and to second 
ones. We will see that the force-constants become less by 
factor 5 while going from the first to the second neighbors 
(see Table 3). Therefore, we do not include the fourth 
neighbors into consideration. 

For the first and third neighbors (in the B sublattice), 
the dynamical matrix has the form 



(q) = 



ELi*.1 B (B K )exp(zqBJ 
=1 ^ B (B' K )exp( iq B;), 



(2) 



them, we introduce variables £,77 = x ± iy transforming 
under the rotation C3 around the z-axis (taken at the A 
atom) as follows (£, 77) — > (£, 77) exp(±27ri/3). In the rota- 
tion, the atoms change their positions Bi — * B2 — + B3, 
Ai — > A 3 — > A 5 , and A 2 — > A 4 — + A 6 . Therefore, 
all the force constants $^ S (B K ) with the different n (as 
well as <f>f B (B K )) are equal to one another, but the force 
constants with the coincident subscripts £ or r\ trans- 
form as covariant variables. For instance, <3?^ B (Bi) = 
$f ( B (B 2 ) exp (27ri/3) = $^ S (B 3 ) exp (-2tu/3). 

The relation between $^ A (A K ) with the points 
Ai,A3,A5 (and also between A4,A2,Ag) has the 
same form. The constants a z — $^(Bi), ^ z = 
a&HAx), < = ^ z B (B[),a = $|f (B x ), a' = 
$^(Bi), and 7 = $£ A (A X ) are evidently real. The 
constant j3 = $^ B (Bi) as well as f3' = ^ B (B[)) is real 
because the reflection (x,y) — > (x, — y) with Bi — * Bi, 
B^ — * B[ belongs to the symmetry group. Besides, we 
have one complex force constant S = $^ A (Ai). 

Two force constants $ AA (A ) and <5> AA (A ) for the 
atom Ao can be excluded in the ordinary way with 
the help of conditions imposed by invariance with re- 
spect to the translations of the layer as a whole in the 
x/z directions. Using the equations of motion JT|) and 
Eqs. ©, ©, we find this stability condition $f r f(A ) + 
6$^ 4 (A 1 )+3$ AB (B 1 ^ • " - ' " 

form for the zz components. 



and the similar 



and for the second neighbors (in the A sublattice) 



Dispersion of bending out-of-plane modes 



^(q) = ^(A )+EL 



)exp(iqAJ,(3) 



where A indices the atom chosen at the center of the 
coordinate system in the A sublattice and the wave vector 
q is taken in units of 1/a. 

The point group D§h of the honeycomb lattice is gen- 
erated by {Cq, a v , er z }, where a z is a reflection z — > — z 
by the plane that contains the graphene layer, C§ is a 
rotation by 7r/3 around the z axis, and a v is a reflec- 
tion by the xz plane. The transformations of the group 
impose constraints on the dynamical matrix. To obtain 



The out-of-plane vibrations u A , uf in the z direction 
are not coupled with the in-plane modes because the force 
constants of type $ IZ or § yz equals zero due to the re- 
flection z — > — z. The corresponding dynamical matrix 
for the out-of-plane modes has the form 



<^(q) <^(q) 

lAB(„\* j,AA 

zz 



'(q)* <^/(q) 



(4) 



where 



A, B 3 




tf> AA (q) = -3(a, + a' z ) 
+27 z [cos(\/3(j a ) + 2cos(3a x /2)cos(V3a y /2) - 3] , 
^ B (q) = a z[exp (iq x ) + 2 exp (~iq x /2) cos (\fSq y /2)\ 
+a' z [exp (-2iq x ) + 2 exp (iq x ) cos (V^q v )] ■ 

(5) 

The phonon dispersion for the out-of-plane modes is 
found 



A 4 B 2 



Wzo,ZA(q) = V^/(q)±|04 B (q)l- 



(6) 



FIG. 1: First, second, and third neighbors in the graphene 
lattice. 



The equations allow us to express the phonon frequencies 
of the out-of-plane branches at the critical points T, K, 
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TABLE I: Force constants in 10 5 cm -2 : a, f3, and a z for the 
first neighbors ; 7, 8, and y z for the second neighbors ; a , a z , 
and /3' for the third neighbors. 



a 


13 


7 


5 


a' 





a z 


7z 


<*; 


-4.095 


-1.645 


-0.209 


0.690 


-0.072 


0.375 


-1.415 


0.171 


0.085 



and M in terms of the force constants: 

^zo(r) = [-6(a 2 +c/ z )] 1/2 
^zo, za(K) = [-3(a, + 4) - 9 7z ] 1/2 
wzo(M) = [-4a, -87,] 1/2 
za(M) = [-2a* - 6a' z - 8 7z ] 1/2 . 



(7) 



Expanding Eq. ([6]) in powers of the wave vector q, 
we find the velocity of the acoustic out-of-plane mode 
propagating in the layer 



, [-0.75a z - 3a' z - 4.5 7z ] 



1/2 



VC^/p, (8) 



where we use the well-known formula for the velocity 
of the acoustic z-mode propagating in the ^-direction in 
terms of the elastic constant C44 and density p of a hexag- 
onal crystal. Because the interaction between the layers 
in graphite is weak, we can correspond the values of C44 
and p to graphite. 



B. Dispersion of in-plane modes 

The dynamical matrix for the in-plane vibrations has 
the form similar to that for the in-plane mode but 
instead of the functions (j)^ (q) and c/>ff (q) we have to 
substitute correspondingly the 2x2 matrices 



TABLE II: Elastic constants (in 10 GPa) and the sound ve- 
locities (in km/s) calculated (theo) and observed (exp) . 





Cn 


C12 


C44 


SLA 


st a 


S z 


theo 


86 


18 


0.57 


19.5 


12.2 


1.59 


exp 


106 ± 2 a 


18±2 a 


0.45 ± .05 a 


w 24 b 


U b 





a Referenced b Referenced^, 



phonon frequencies for the in-plane branches at T and K 
are found 



u> 



in— pi 1 
1,2 



(r) = [-6( a + a')] 1/2 , 

..|7 pl (^) = [-3(a + a')-97] 1/2 , 
(K) = [-3(q + a') - 9 7 ± 3(/3 - 



in— pi 1 
jn-pl, 



doublet, 
doublet, (10) 
-P')] 1/2 - 



An algebraic equation of the forth order have to be solved 
in order to find the phonon frequencies at the M point as 
well as at points of the general position. 

The in-plane vibrations make a contribution into the 
elastic constants Cn and C12. The corresponding rela- 
tion between the dynamic matrix elements and the elas- 
tic constants can be deduced taking the long-wavelength 
limit (q — > 0) in the matrices ©■ In this limit, separat- 
ing the acoustic vibrations u ac from the optical modes, 
we obtain the equation of motion in the matrix form 



^AA + (j) AB + (j) BB + (j) BA )/2 



0. 



(ii) 



where the subscripts and 1 mean that the terms of the 
zero and first order in q should, correspondingly, be kept 
in the matrices ([9]) , but the expansion to the second order 
is used in other terms. We find the matrix factor of u ac 



^(q) </>f/(q) 
<^(q)* ^(q) 



^f(q) ^f(q) 



<C(q) <Oi) 



(9) 



The matrix elements ^^(q) and <£^(q) are obtained 
from ^p(q) and (q), Eqs. (0, correspondingly, 
with substitutions 7, a, and a' instead of j z , a z , and 
a' z . The off-diagonal elements are given by 

^ A (q) = 

S[exp(iy3q y ) + 2 cos (3q x /2 + 2tt/3) exp(— iy3q y /2)] + 
S* [exp(-iy/3q y ) + 2 cos (3q x /2 - 2tt/3) exp (iVSq y /2)], 

<C(q) = 

/3[exp (iq x ) + 2 exp (-iq x /2) cos (y/3q y /2 - 2n/3)] 
+/3'[exp (—2iq x ) + 2 exp (iq x ) cos (y/3q y + 27r/3)]. 

The matrix elements for the B sublattice can be obtained 
from that for the A sublattice by C2 rotation (x,y) — > 
— (x,y) of the graphene symmetry group. The optical 




M K r 

reduced wave vector 

FIG. 2: Calculated phonon dispersion for graphene; the force 
constants, elastic constants, and phonon frequencies at critical 
points are listed in Tables 1, 2, and 3 correspondingly. 
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TABLE III: Phonon frequencies at critical points in cm ; z and |j stand for the out-of-plane and in-plane branches, respectively. 





r 


[0 0] 


M [1 V3]7r/3a 


K [0 l]47r/3V3a 
















wf 


II II II 

, ,11 , ,11 . ,11 , .2 

U l ^2,3 ^4 ^1,2 


exp 


loyu 


oDl 


1389" 




630 d 


670" 


471 c 


lolo _Llo4 4oZ 




1583 6 


868 c 


1390 6 


1323 ft 1290 6 






451 d 


1265 ft 1194 6 517 d 




1565 6 


868 e 






625 e 


625 e 


480 e 


1285 e 1021 e 537 e 


theo^ 


1595 


890 


1442 


1380 1339 


636 


618 


475 


1371 1246 994 535 


theo 6 


1581 




1425 


1350 1315 








1300 1220 


theo 


1581 


893 


1363 


1324 1279 


651 


655 


308 


1349 1199 1028 495 



Reference^, b Reference^, c Reference^, d Reference^, e Reference^, ^ Reference^ 



in Eqs. (fTTj) : 

( Siq2 - Uj2 tO, (12) 

\ S2?i sir - ^ ) 

where 

si = -§7 - fa - 3a' + - 2f3') 2 /{a + a'), 
s 2 = f Re(5) - |/3 - 

With the help of Eq. (| 1 2[) . we obtain the velocities of 
longitudinal and transverse acoustic in-plane modes 

s L A = a^Jsi + s 2 = a/Ch//j, 
s T a = aV s i _ s 2 = yJflxx ~ Ci 2 )/2p, 

corresponding them to the elastic constants Cn, C12 and 
density p of graphite. 

III. RESULTS AND DISCUSSIONS 

The calculated phonon dispersion is shown in Fig. 2. 
Notice, first, that the sound velocities (for long waves, 
q — * r) are isotropic in the xy plane as it should be 
appropriate for the symmetry of graphene. Second, the 
in-plane LO/TO modes at T, the in-plane LO/LA modes 
at K, and the out-of-plane ZA/ZO modes at K are dou- 
bly degenerate because graphene is the non-polar crystal 
and the C^ v symmetry of these points in the Brillouin 
zone admits the two-fold representation (observation of 
splitting of that modes in graphene would display the 
symmetry braking of the crystal) . 

Because of the lack of information on graphene, 
we compare the present theory with experiments on 
graphite. Thus, we have only three force constants a z , 
j z , and a' z to fit four frequencies of the out-of-plane 
modes at the critical points T, M, and K. We must 
keep in mind that the frequencies in graphene for the 
out-of-plane branches could be less than their values in 
graphite since the atoms are more free to move in the 
z direction in graphene comparatively with graphite. It 
is evident that the adjacent layers in graphite affect the 
low frequencies more intensively. The interaction of the 



adjacent layers can be estimated from the ZA - ZO' split- 
ting about 130 cm -1 given, for instance, in Ref^. These 
modes become degenerate when the inter-layer interac- 
tion is switched off. Therefore, the lowest frequencies of 
out-of-plane modes calculated at the M and K points 
are considerably less than the corresponding frequencies 
observed in graphite (see Table 3). 

Furthermore, the force constants determine the veloc- 
ity s z , Eq. ([5]), of the acoustic out-of-plane mode along 
with the elastic constant C44. We see that the velocity 
has the nonzero value unless a definite condition is sat- 
isfied for the force constants. Using the values of force 
constants obtained in fitting with the experimental data 
(see Table 1), we find the value of the sound velocity 
s z = 1.59 km/s for the out-of-plane mode. This result is 
contradictory to the statement of ReffSS, that the acoustic 
out-of-plane mode has a quadratic dispersion. The fact 
that the sound velocity s z is very sensitive to the small 
variation of 7 Z indicates that graphene is nearly unstable 
with respect to transformation into a phase of the lower 
symmetry group at T. 

For the in-plane modes, we have to fit eight frequencies 
at the critical points and two elastic constants. Equations 
(dU]) and (fl"3"|) can be used as a starting point. Fitting of 
the in-plane branches is insensitive to the imaginary part 
of the constant S. Therefore, it is taken as a real parame- 
ter. Results of the fit are presented in Fig. 2 and Tables. 
Notice, that the extent of agreement of the present theory 
with the data obtained for graphite corresponds to the 
comparison level between the first-principle calculations 
for graphite in Ref^ and their experimental data (see 
Table 3). The largest disagreement of 5% between our 
calculations and experiments on graphite for the highest 
phonon mode occurs at the K point. This is result of the 
Kohn anomaly due to the electron-phonon interactional 
which reduces the phonon frequency at K. The same 
reason explains some overbending observed probably in 
graphite along the T — M direction. 

IV. CONCLUSIONS 

We calculate the phonon dispersion in graphene us- 
ing the Born- von Karman model with the first-, second-, 
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and third-neighbor interactions imposed by the symme- 
try constraints. The bending (out-of-plane) modes are 
not coupled with the in-plane branches and indicate the 
latent instability of graphene with respect to transforma- 
tion into a lower-symmetry phase. The acoustic ZA mode 
has the linear dispersion in a small wave-vector interval 
near the T point. The optical frequencies of these modes 
are less than the corresponding values in graphite. For 
the higher in-plane modes, the fit shows good agreement 



between the experimental and calculated values of optical 
frequencies, elastic constants, and acoustic velocities. 
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